Collisions in zero temperature Fermi gases 
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We examine the collisional behavior of two-component Fermi gases released at zero temperature 
from a harmonic trap. Using a phase-space formalism to calculate the collision rate during expan- 
sion, we find that Pauli blocking plays only a minor role for momentum changing collisions. As a 
result, for a large scattering cross-section, Pauli blocking will not prevent the gas from entering the 
collisionally hydrodynamic regime. In contrast to the bosonic case, hydrodynamic expansion at very 
low temperatures is therefore not evidence for fermionic superfluidity. 
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The last few years have seen rapid progress in the field 
of ultracold atomic Fermi gases [1-6]. Most recently, 
regimes of strong interactions have been observed in these 
gases near Feshbach resonances [7-10]. Studies of these 
systems are of particular importance because of the pos- 
sibility of creating BCS-likc supcrfiuids [11,12]. Such a 
realization would establish highly controllablc model sys- 
tems for studying novel regimes of fermionic superfluid- 
ity 

A unique feature of atomic systems is the ability to 
analyze the gas by turning off the trapping potential 
and observing the expansion. The expansion behavior 
can reveal the momentum distribution and the effects 
of mean- field interactions and collisions. Hydrodynamic 
behavior can be easily detected when the gas is released 
from an anisotropie atom trap. In that case, the spatial 
anisotropy of the cloud reverses during free expansion. 
This is caused by the larger pressure gradient along the 
tightly confining direction, which leads to a faster expan- 
sion, and subsequent reversai of the spatial anisotropy. 
This anisotropie expansion was used to identify the for- 
mat ion of the Bose-Einstein condensate [13,14]. 

A BEC obeys the hydrodynamic equations of a su- 
pcrfluid [15]. However collisional hydrodynamics arising 
from a high elastic collision rate also results in anisotropie 
expansion [16,17]. Thus, the normal component can also 
expand anisotropically [18]. For the bosonic case, two key 
points make the distinction between the two fractions ob- 
vious: (i) At the typical transition temperature, the BEC 
has much less energy than the normal cloud, so the two 
components are clearly separated in size. (ii) The scat- 
tering rate needed to achieve condensation is usually not 
large enough that the normal gas is in the hydrodynamic 
regime. For these two reasons, the appearance of a dense 
anisotropie cloud during expansion is considered to be 
the "smoking-gun" for the formation of a Bose-Einstein 
condensate. 

A superfluid Fermi gas is predicted to obey the su- 
perfluid hydrodynamic equations of motion [19-22] and 
therefore should show strong anisotropie expansion when 
released from an anisotropie harmonic trap [22]. The re- 
cent observation of anisotropie expansion of an ultracold, 



interacting, two-spin fermionic mixture [7-9] has created 
considerable excitement and raised the question under 
what conditions is this expansion a signature of fermionic 
superfluidity and not of collisional hydrodynamics. There 
are two major differences from the bosonic case: (i) Since 
the energy of ultracold fermions always remains on the 
order of the Fermi energy, the size in expansion for both 
normal and superfluid components will be similar. (ii) 
Current efforts towards inducing BCS pairing ali take 
place in strongly interacting systems. This results in a 
large scattering rate modified only by the effects of Pauli 
blocking at low temperatures. 

The interpretation of the observed anisotropie expan- 
sion in strongly interacting Fermi gases is therefore crit- 
ically dependent on the role of Pauli blocking of colli- 
sions during the expansion. The tentative interpretation 
of anisotropie expansion as superfluid hydrodynamics [7] 
was based on the assessment that collisions are strongly 
suppressed at sufficiently low temperatures [7,23-28]. 

Here we show generally that the collision rate becomes 
independent of temperature and prevails even at zero 
temperature, if the Fermi surface is strongly deformed. 
This happens in an extreme way during ballistic expan- 
sion. In the small cross-section limit, we find that less 
than half of the total number of momentum changing col- 
lisions is suppressed. For a large scattering cross-section, 
the absence of suppression results in strong collisional 
behavior of normal Fermi gases during expansion for ali 
initial temperatures. This result has the important con- 
sequence of rendering expansion measurements of Fermi 
gases near Feshbach resonances ambiguous for differenti- 
ating between superfluid and normal components. 

We first consider the expansion of a single component 
Fermi gas. At ultralow temperatures, fermionic anti- 
symmetry prevents s-wave scattering in a single compo- 
nent and renders the gas completely collisionless. The 
phase space occupation f{x\, X2, x$,Pi,P2,P?ì) = /(x,p) 
at zero temperature in a harmonic trap with frequencies 
CJ2, W3) can be written as 

/(x,p) = Q{E F - Eimu?a$/2 - Z lP 2 l /2m) 

where m is the particle mass, is the Heaviside step 
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function defìned as 9(x) = 0(1) for x < 0(x > 0) and 
E F = hi&NIi.uJi) 1 ^ is the Fermi energy for N particles. 
At time t = 0, the trapping potential is turned off sud- 
denly, allowing the gas to expand freely. At t — 0, the 
momentum space Fermi surface at x = (xi, x%, x%) is 



Y>ipj/2m = Ep — T,iTnw i x i /2, 



(1) 



a sphere of radius \j2ra{Ep — 'Siinujfxf/2). In this non- 
interacting system, the evolution of the Fermi surface 
can be derived from the simple evolution law for ballistic 
expansion x(0) = x(t) — pt/m. Substituting this in Eq. 1, 
we obtain: 



local-density approximation. These collisions have an ap- 
pealing geometrical picture in the locai phase-space de- 
scription (Fig. l(c)). Each elastic collision involves one 
particle from each spin state. We label with p's and q's 
the momenta of the two different spin states. Considcr 
the collision Pi + qi — > P2 + q2- Conservation of momen- 
tum and kinetic energy mandates P2 + q2 = Pi + qi and 
|P2 — qa| = |pi — qi |. These relations restrict P2 and q2 
to lie on diametrically opposite ends of the sphere with 
Pi — qi as a diameter. The deformation of the Fermi sur- 
face during expansion opens up unoccupied final states 
P2 , q2 and therefore allows collisions to take place even 
in a zero temperature Fermi gas (Fig. l(c)). 



2m 



uift 



2,2 



t 1 



H 2 



2 2 

mujtxt 



2(1 



(2) 



which is an ellipsoid with generally unequal axes 
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The anisotropy of the Fermi surface during expansion can 
be understood generally by noting that for long times t, 
at any position x, the spread in momentum Api(t) can 
only arise from the initial spread in position Axì(0). For 
anisotropie traps this gives rise to an anisotropie mo- 
mentum distribution during ballistic expansion. For a 
mixture of two spin states, this deformation of the Fermi 
surface from a sphere into an anisotropie ellipsoid re- 
moves Pauli blocking of final states and allows collisions, 
as will be shown. 

The momentum distribution at position x given by 
Eq. 2, also allows us to simply calculate the spatial den- 
sity distribution as the volume of the momentum-space 
ellipsoid, 



2mE F \ / i 1 ~ 2E^^i i+Jh* l 



h- 



1^(1+^2)1/2 



(3) 



in agreement with other derivations [29] . For long expan- 
sion times t, the spatial distribution becomes isotropie, 
mirroring the isotropie momentum distribution in the 
trap. 

Specializing to the experimentally relevant case of 
a cylindrically symmetric trap, ballistic expansion de- 
forms the locai Fermi surface into a momentum ellip- 

soid of cylindrical symmetry with aspect ratio w 1+UJ ì p 

(Fig. l(a,b)). Here oj±(uj z ) is the radiai (axial) trapping 
frequency. For long times t, this deformation reaches the 
asymptotic aspect ratio u> z /u>± = A, the initial spatial 
aspect ratio in the trap. 

Now consider an equal mixture of two spin states which 
interact via a finite s-wave scattering length. We assume 
that the trapping frequencies are identical for the two 
states (standard experimental conditions) and special- 
ize to the usuai case of two-body elastic collisions in the 
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FIG. 1. (a) Deformation of the momentum space Fermi sur- 
face at x = 0, from a sphere to an ellipsoid during expansion 
from an anisotropie harmonic trap. The case of cylindrical 
symmetry is shown, where the three-dimensional distribution 
is symmetric about the vertical axis. The parameters chosen 
are an aspect ratio A = 0.2 and expansion times u>±t = 0,1,3 
and 10. (b) The deformation at a position radially displaced 
by Ef /mu) 2 ± . (c) Geometrical representation of collisions 
in momentum space. The two spin states have identical dis- 
tributions. Three different types of collisions are shown for 
particles with initial momenta pi and qi - none, one or both 
of the final states are occupied. 

The effect of collisions can be formally calculated from 
the Boltzmann transport equation for the evolution of 
the phase space distribution /(x, p, t). In the absence 
of external potentials and ncglecting mean field we have 
[30]: 



df df 
at ox 



(4) 



where v = p/m and r co ii[/] describes the effect of col- 
lisions. Collisions attempt to restore locai equilibrium 
by countering the deformation of the momentum space 
Fermi surface during free expansion (Eq. 2). 
r c oii[/] can be written as the collision integrai: 



r(x, Pl ,i) = - i ^ r ; d 3 qi d*n 



|pi-qi| 



[/(pi)/(qi)(l-/(P2))(l-/(q 2 )) 
-/(p 2 )/(q 2 )(l-/(pi))(l-/(q 1 ))] 



(5) 
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where a is the momentum-independent scattering cross- 
section, /(pi) = /(x,pi,t), /(q ; ) = /(x,q ; ,i) and fi 
points along p 2 — q_2 ■ The integrai over qi is over the mo- 
mentum ellipsoid at position x and time t for one of the 
spin states. The first term in the integrand is the collision 
rate for the process Pi +qi — * P2 + q2 The second term 
corresponds to the reverse process P2+q2 — ► Pi+qi, and 
ensures that only distribution changing collisions con- 
tributo . 

Pauli blocking is expressed in the suppression factors 
for the final states (1 — /) in V. The collision integrai 
neglecting Pauli blocking, Tci.p, is furnished by setting 
the suppression factors ali equal to 1 in Eq. 5. This is 
the rate for classical collisions which change the momen- 
tum distribution. The total classical collision rate Tei 
is the first term on the right hand side of Eq. 5 without 
any suppression factors. In addition to Tci.p, this also 
contains the rate for collisions which do not change the 
momentum distribution: if both final states are occu- 
pied, then the reverse process has the same rate. These 
additional collisions do not affect observables of the sys- 
tem. Fig. l(c) shows examples of these different types of 
collisions. Pi + qi — > P2 + q2 contributes to Tei, Tchp 
and r. Pi + qi — > P3 + q3 contributes to Tei and rci, p . 
Pi + qi — > P4 + q4 contributes only to Tei- To deter- 
mine the effect of Pauli blocking, we compare Y and Tci.p 
for a small cross-section a. The collision rate at a par- 
ticular time t can then be calculated perturbatively, by 
propagating the system ballistically for the time t and 
then evaluating Eq. 5 with and without the suppression 
factors. 

Fig. 2(a) displays the numerically calculated collision 
rates T, rci, p and Tei, evaluated at x = p = 0, as a 
representative case, for an initial aspect ratio A = 0.03. 
Both r and rci, P increase initially as the deformation 
of the Fermi surface becomes more pronounced. For 
long times (uj±t ^> 1), they are both suppressed be- 
cause both the density (J d 3 qi) and the relative ve- 
locity ( ^ Pl m qi ) drop. The two curves approach each 
other with time since Pauli blocking becomes less ef- 
fective with stronger deformation. The fraction of mo- 
mentum changing collisions which are not affected by 
Pauli blocking, F(X) = J dtT(0,0,t)/ J dtT C i, P (0,0,t), 
is shown in Fig. 2(b). The main result of our paper is 
the inefficiency of Pauli blocking during expansion from 
anisotropie traps. For A < 0.05, F > 0.5, and approaches 
~ 0.55 as A approaches 0. Most experiments work in this 
regime of trap aspect ratio. 

The above results form an upper bound on the Fermi 
suppression even if we consider ali the possible collisions 
occurring in the system, for arbitrary x and p. First, we 
observe that for ali x, at any time t, the Fermi surface 
is identically deformed and different only in size accord- 
ing to the locai density (Eqs. 2,3, Fig. l(a,b)). We have 
checked numerically that to within 5%, the centrai mo- 



mentum provides a lower bound on T within a momen- 
tum ellipsoid, at ali x and for ali times t. Next, we note 
that for x ^ 0, the density n(x, t) puts more weight at 
longer times than n(0, t) (Eq. 3). Since Fermi suppres- 
sion becomes less cffcctivc with time, Pauli blocking is 
most cffcctivc at x = 0. The calculation for x = p = 
thus provides an effective upper bound for the overall 
collisional suppression in the system. We conclude that 
more than half of ali the possible collisions are not Pauli 
blocked for typical experimental values of A. 
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FIG. 2. (a) Collision rate as a function of expansion time 
in the perturbative approximation for the initial aspect ratio 
A = 0.03. Dashed line - total classical collision rate r c i, thin 
line - classical rate for momentum changing collisions r c i lP , 
thick line - collision rate for fermions F. The displayed rates 
were evaluated at x = and p = and give an effective upper 
bound on the Fermi suppression. (b) Allowed fraction of col- 
lisions F(X) for a zero-temperature two-spin Fermi gas. For 
an initial aspect ratio A = 0.05, F is 0.5. For large anisotropy 
(A — » 0), F approaches ~ 0.55. 

So far, we have not considered the effect of the colli- 
sions themselves on the momentum distribution. Colli- 
sions drive the system towards equilibrium, which cor- 
responds to an isotropie Fermi-Dirac distribution. If 
this collisional relaxation (Eq. 5) is much faster than the 
non-equilibrium perturbation due to ballistic expansion 
(Eq. 2), the momentum distribution maintains locai equi- 
librium at ali times. If locai equilibrium is maintained, 
the Boltzmann equation leads to the hydrodynamic equa- 
tions [30]. For free expansion from anisotropie atom 
traps, these equations lead to the reversai of anisotropy 
[16,17]. Even if equilibrium is not fully maintained, col- 
lisions always have the effect of transferring momentum 
from the weakly confining axis to the strongly confining 
axis resulting in an eventual spatial aspect ratio > 1 [17]. 

We now want to reconcile our new result that Pauli 
blocking is inefficient during free expansion, with pre- 
vious results [26] which show that at low temperatures, 
collisional damping of collective excitations is suppressed. 
For this, we derive an equation of motion for the momen- 
tum space anisotropy a to leading order in a and T/Tp 
[31]: 

à = \{d x v x + d y Vy - 2d z v z ) - — —C (a, — ] (6) 
o m \ 1 p I 
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where C describes the collisional relaxation and has the 
asymptotic forms: 
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a »{T- F )- ™ 
In terms of a, the aspect ratio of the momentum space 
ellipsoid is yj l+2a ■ PF,T,Tp are the locai Fermi momen- 
tum, temperature and Fermi temperature respectively. 
This equation was derived from the second momentum 
moment of the Boltzmann equation (Eqs.4,5), using a 
Fermi-Dirac distribution with an anisotropie Fermi sur- 
face as ansatz [32] . The numerical coefhcicnts in Eq. 7 
were obtained by analytic integrations over momentum 
space. 

At zero temperature, there is no linear term in a in 
Eq. 7. This shows that Pauli blocking is efficient as long 
as the anisotropy is small. This is the case for small am- 
plitude excitations in a trapped degenerate gas [26] . How- 
ever, for the large anisotropies of ballistic expansion, the 
a 3 term, which is independent of temperature and not 
affected by Pauli blocking, is responsiblc for collisional 
relaxation. 

Eqs. 6,7 allow us to distinguish collisionlcss from hy- 
drodynamic behavior in different regimes. The driving 
term involving v is on the order of the trap frequency 
uj± and the damping term has a prefactor navp. Thcrc- 
forc, the dimensionlcss parameter characterizing the at- 
tainment of the hydrodynamic limit is <&o = navp/uJ±. If 
$o <C 1, then one can neglect collisions entirely, and the 
gas will expand ballistically. For small anisotropies, hy- 
drodynamic behavior requires <&o(T/Tf) 2 3> 1. For large 

1/3 

anisotropies, hydrodynamic behavior requires $ 3> 1. 
At ultralow tempcratures, the expansion after release 
from a highly anisotropie trap may be collisionlcss ini- 
tially, but as a grows, the a 3 term in Eq. 7 will become 
important, and induce hydrodynamic behavior. 

Our calculations clearly predict that for parameters 
of current experiments, <&o > 1, free expansion will not 
be collisionless, but show behavior which is at least in- 
termediate between collisionless and hydrodynamic [17]. 
Full hydrodynamic behavior may not be achieved, since 
for small values of a, Pauli suppression becomes effective 
again. More quantitative studies are necessary in order 
to assess how much this behavior would differ from su- 
pcrfluid expansion. This could be realized by extending 
analytical studies [17] to high degeneracies or by Monte- 
Carlo techniques [33]. Our main conclusion is clcar, how- 
evcr, that the breakdown of Pauli blocking under free 
expansion means that hydrodynamic expansion will not 
be the dramatic, qualitative signal for superfluidity in 
strongly interacting fermions, the way it was for BEO 
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